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ABSTRACT 

We present the results of an analysis of charged particle pseudorapidity 
distributions in the central region in collisions of a Au projectile with Al, 
Cu, Au, and U targets at an incident energy of 10.8 GeV/c per nucleon. 
The pseudorapidity distributions are presented as a function of transverse 
energy produced in the target or central pseudorapidity regions. The corre- 
lation between charged multiplicity and transverse energy measured in the 
central region, as well as the target and projectile regions is also presented. 
We give results for transverse energy per charged particle as a function of 
pseudorapidity and centrality. 

PACS number(s): 13.85.-t, 25.75. +r 

1 Introduction 

Global observables such as (pseudo-) rapidity particle density contain valu- 
able information on the reaction dynamics and, indirectly, on the degree of 
thermalization as well as the energy and entropy densities reached in rel- 
ativistic nucleus-nucleus collisions. With light projectiles, pseudorapidity 
distributions have been studied in detail previously both at AGS and CERN 
energies (for a review, see Ref. [1]). Our collaboration earlier reported the 
measurements of charged particle distributions , energy flow and stopping 
[Q], transverse energy distributions |Q, in the collisions of a Si beam with 
Al, Cu, and Pb targets. Large energy deposition has been inferred from 
these measurements. Extrapolating these results to heavy projectiles raises 
expectations to create, in these collisions, the deconfined phase of quarks 
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and gluons. We report here measurements of the charged particle pseudo- 
rapidity distributions in colhsions of 10.8- A GeV/c Au beams with several 
nuclear targets carried out in Experiment 877 operating at the Brookhaven 
National Laboratory Alternating Gradient Synchrotron (AGS). We combine 
the multiplicity data with our earlier measurements 0] of transverse energy 
(Et) pseudorapidity distributions in Au+Au collisions at the same energy to 
study the Et per charged particle as a function of pseudorapidity and the 
centrality of the collision. 

2 Experimental setup 

The E877 experimental setup is shown schematically in Fig. |l[ For the pseu- 
dorapidity distribution analysis we use data primarily from the Multiplicity 
Detector, complemented by data from the Participant and Target Calorime- 
ters (see the insert in Fig. |1|), which provide a measurement of the centrality 
of a collision. "Zero-degree" energy (deposited mainly by projectile spec- 
tators) is measured by the Uranium Calorimeter, situated in the forward 
spectrometer. The horizontal position of the beam particle is measured by 
a pair of silicon strip detectors, the Beam Vertex Detectors, shown in the 
insert in Fig. |I| (BVer 1 and BVer 2). The information from these detectors 
is used on an event by event basis. The mean vertical displacement of the in- 
coming beam particle is estimated using the information on the distribution 
of the "centroid of hits" in the multiplicity counter (see details in Section 
1^). The angular divergence of the beam (~ 1 mr) is much smaller than the 
bin widths in r] and (p used for multiplicity distribution analysis. Most of 
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upstream interactions are effectively rejected using the pulse height from a 
100 fim thick Si surface barrier detector located just upstream of the target. 

The Multiplicity Detector, shown in Fig. ^ consists of two identical silicon 
pad detectors, each of which was made from a disc of silicon 300 fim thick 
and approximately 3.8 cm in radius. To reduce the number of 5-electrons 
reaching these detectors, two 3 mm thick aluminum absorber plates were 
placed upstream of each plane of silicon. The active region of each detector 
is a ring of inner radius 1.4 cm and outer radius 3.4 cm divided into 512 
pads. The detectors are segmented into 8 concentric rings of 64 pads each. 
One detector, located 3.37 cm from the target, covers the pseudorapidity 
region 0.87 < r] < 1.61, and the other, located 8.17 cm from the target, 
covers the region 1.61 < rj < 2.46. These values of pseudorapidity coverage 
correspond to the case where the beam particle is incident at the center of 
the detector. Due to the finite size of the beam spot and variation of the 
beam position during the AGS spill, the actual pseudorapidity coverage is 
slightly larger. The size of the pads corresponds approximately to 0.1 in 
both 7] and azimuthal angle 0, which determines the angular resolution in 
this measurement. Signals from the pads, after preamplification and shaping, 
are sampled at the peak and digitized. For the most central events the mean 
occupancy in pads which see the highest track density is close to 0.3. 

The Participant Calorimeter (PCal) ||^ is a lead/iron/scintillator sam- 
pling calorimeter. It has a depth of four interaction lengths and a radius 
of approximately 84 cm. It is approximately azimuthally symmetric, built 
with four identical quadrants. Each quadrant of the PCal is divided into 
four azimuthal slices of 22.5". Each slice is divided radially into eight towers. 
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Longitudinally, the calorimeter is divided into two electromagnetic depth 
segments and two hadronic depth segments. This division leads to a to- 
tal of 16x8 = 128 towers for each quadrant and 512 towers for the entire 
calorimeter. PCal measures energy flow into the polar angle region which 
corresponds to pseudorapidity range 0.83 < rj < 4.7. The Target Calorime- 
ter (TCal) is made of 992 Nal crystals each 5.3 radiation length deep. It 
covers the backward hemisphere, corresponding to the pseudorapidity range 
—0.5 < f] < 0.8. For more details on TCal and the analysis of TCal data 
see The Uranium Calorimeter (UCal) consists of 25 modules and mea- 
sures the energy of particles entering the forward spectrometer through a 
collimator with an opening of -115 mr< 6^ <14 mr and -21 mr< 6y <21 mr. 

The data were taken with several targets, Al (242 mg/cm^, approximately 
1.9% of an interaction length for a gold projectile), Cu (500 mg/cm^, ~2%), 
Au (540 mg/cm^ and 980 mg/cm^, ^1% and 1.8%), and U (575 mg/cm^ and 
1150 mg/cm2, ^1% and 2%). 

3 Analysis 

3.1 Pulse height spectra 

The pulse height distributions in the Multiplicity Detector were first cor- 
rected for pedestal offsets and differences in gain, and the non-functional 
channels were identified. A channel was defined as good if the pulse height 
distribution showed a minimum ionizing particle (m.i.p.) peak well separated 
from the pedestal (as in the example shown in Fig. Dead or noisy chan- 
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nels (altogether about 25%) were removed in the analysis. The corrected 
pulse height distributions were studied in a variety of ways. In the vicinity of 
pedestals the distributions were fitted by a Gaussian; this fit gives the width 
of the electronic noise distribution in each particular channel, and permits 
to evaluate the mean occupancy of the pad (1 —po, where po is the probabil- 
ity of the pad not being occupied). The part of the distributions above the 
pedestals were fitted by a Landau distribution convoluted with a Gaussian 
describing the electronic noise and taking into account the effects of atomic 
binding of the electrons P| (this fit is similar to the analysis done in 0): 

/(A) = r hie) exp( ~^^ 7^' )& (1) 

^ = (52 + aLJ^/^; (2) 

Here A is the actual energy loss, the variance 62 is related to the effect of 
electron atomic binding, and a^gj^g is the variance of electronic noise; fii^) 
is the Landau distribution function: 

1 1 rc+ioo 

hie) = -$(A); <I>(A) = — / exp(n + \nu + \u)du (3) 

A = |(e - (emp - eAo)) = '-^^ + Ao; Aq = -0.225, (4) 

where ^ is the width of the distribution, and e^p is the most probable energy 
loss; c is an arbitrary real positive constant. For the case of multiple (n) hits 
the parameters ^ and e^p are to be replaced by: 

a = < (5) 
en,mp = n{e^p + ^ ln{n)) (6) 
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To investigate the distribution in number of hits for the purpose of eval- 
uation of mean pad multiphcity, two different fits to the pulse hight distri- 
butions were carried out. The first fit assumed a Poisson distribution in the 
number of hits in the pad, and the parameter extracted was the mean occu- 
pancy, which can be compared with the value extracted from the fit to the 
electronic noise distribution. The quality of the fit can be seen in Fig. ^ 
where the fitted curves are shown along with the data from one of the pads. 

In the second fit, the probabilities of single or double hits are free pa- 
rameters. We find that the probability of double hits defined independently 
exceeded the value expected from Poisson statistics. For pads with mean 
occupancy ~ 0.3 the ratio of observed double hits to the calculated value 
from Poisson statistics is ~ 1.3. This effect is understood to be the effect of 
7-conversions in the target and in the absorber. The Monte Carlo simulation 
(described below) shows that for a heavy target (Au or U) about 5% of all 
hits are due to 7's from 7r° decays converting into e^e^-pairs. Due to the 
small opening angle about 30% of produced e"'"e~-pairs occupy the same pad 
of the multiplicity detector, introducing a non-Poissonian element in the dis- 
tribution. The correction for 7-conversions used in the analysis which takes 
this effect into account is described below. 

After the calibration was done and the procedure for the calculation of 
the mean pad multiplicity was established, a "hit" threshold was introduced 
corresponding approximately to one-half of the peak of the minimum ionizing 
particle signal. By varying the threshold position and comparing the results 
with the information from the fit to the electronic noise distribution, it was 
found that the occupancy can be defined in this way with an accuracy better 

7 



than 2%. 



3.2 Corrections and Selection Criteria 

Event selection. In order for an event to contribute to the multiphcity anal- 
ysis, it is necessary to obtain the horizontal position of the beam particle. 
Thus events with missing or ambiguous information from the vertex detector 
were rejected. Also, it is important to reduce the background from interac- 
tions upstream of the target as much as possible. Therefore, wc reject events 
in which the pulse height in the upstream silicon detector is below a threshold 
value close to the energy loss peak of Au ions. 

Beam position. To calculate the pseudorapidity corresponding to each 
pad the knowledge of the position of the interaction point relative to the 
multiplicity detector is very important. As mentioned above, the horizontal 
position is measured for each event in the Beam Vertex Detector. The fol- 
lowing technique was used to define the relative position of the Beam Vertex 
Detector with respect to the Multiplicity Detector. We exploit the fact that 
the hit centroid distribution is expected to be axially symmetric when the 
beam particle is incident at the center of the Multiplicity Detector. In each 
event the horizontal and vertical components of the hit centroid position are 
given by 

Ch^YI cos 0i; Cv ^Yl sin (7) 

i i 

where the sum is over all hits and 0j is the azimuthal angle of the pad 
containing the i-th hit. To avoid a bias in this part of the analysis from dead 
pads, a symmetrized dead pad mask was used, declaring some good pads as 
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dead to make the distribution of dead pads symmetric. After the coordinates 
of the centroid are calculated, the horizontal component is plotted against 
the position provided by the Beam Vertex Detector. The results of this 
analysis are shown in Fig. ^. The beam position from BVer corresponding to 
Ch = gives the the relative BVer and multiplicity detector displacement 
(about 0.75 mm from Fig. |). To understand the statistical fluctuations in 
the plot near x = 0, note that the most probable beam position from the 
BVer is about —4 mm; thus the beam position near zero is very rare. 

The average value of Cy gives information on the mean vertical position of 
the beam for any event sample. A plot of Ch versus the horizontal coordinate 
obtained from BVer fixes the scale factor between the hit centroid position 
and the beam displacement. For the data shown in Fig. ^ this scale factor 
is 4.1 mm, which gives us about 1.2 mm for the mean vertical beam offset. 
One can also see from the plot that there is no correlation between vertical 
and horizontal components of the beam position. 

The correction for multiple hits was done pad by pad using the average 
occupancy. Since the mean pad occupancy depends on beam position, the 
correction was done separately for each value of this variable. For this pur- 
pose the corresponding pseudorapidity, azimuthal angle, and solid angle in 
these variables were calculated for every pad as a function of both vertical 
and horizontal beam position. The Multiplicity Detector consists of two sili- 
con pad detectors located at different distances from the target and covering 
different pseudorapidity regions. For non-zero beam positions these regions 
overlap. It was verified that the distributions calculated using the data from 
different detectors coincide in the overlapping region for the values of the 
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beam mean vertical displacement and the Multiplicity Detector relative off- 
set used in the analysis. The other independent check is that the resulting 
azimuthal angular distributions are ffat in different pseudorapidity windows. 
These methods are sensitive to displacements at the level of 0.2-0.3 mm. 

Upstream interactions. Although most of the upstream interaction events 
are effectively rejected using the information from the upstream Si detector 
and the Beam Vertex Detectors, it is important for collisions of medium 
centrality and for data taken with light targets to perform a background 
subtraction to correct for residual upstream interactions. The subtraction 
was done using the data taken with empty target frame. The relative con- 
tamination by upstream interactions is different for events triggered by TCal 
or PCal. Almost all upstream interactions result in relatively low TCal Et, 
but PCal Et can be rather significant. For heavy targets the admixture 
of upstream interactions in the event sample is negligible for events with 
TCal Et > 10 GeV, and only about 20% for the region Et ^4-5 GeV (the 
corresponding differential cross section da/dEt is presented in the following 
section in Fig. H). For the Al target the admixture is about 40% in the same 
region (TCal Et ~4-5 GeV). In contrast to the heavy targets the admixture 
of upstream interactions for the Al target increases with TCal Et because of 
the sharp drop of da/dEt for interactions in the target; for events with TCal 
Et > 7 GeV, upstream interactions become dominant. For the Au target the 
admixture of upstream interactions in PCal Et regions centered at 65 GeV, 
130 GeV, and 190 GeV is about 70%, 40%, and 10% respectively. 

Delta electrons. One of the most important corrections is the subtraction 
of hits due to 5-electrons produced in the target. The aluminum absorber 
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located in front of each detector plane reduces the number of 5-electrons by 
about a factor of 10, but nevertheless their contribution is not negligible. It is 
not possible to extract the pseudorapidity distributions of 5-electrons (needed 
for the correction) directly from the data. The data for "beam" events (no 
interaction in the target) provide the number of produced 5-electrons and 
their pseudorapidity distribution for the case when the incoming nucleus 
traverses the entire length of the target. Unfortunately, this distribution 
cannot be used for the correction for normal events, because the path length 
of the incoming nucleus in the target is different from that of a beam track, 
and thus the multiple scattering and absorption effects are also different. 
To understand the effect of the absorbers on the (5-ray energy and angular 
distribution and to calculate the pseudorapidity distribution needed for the 
correction a detailed simulation was performed using the GEANT (version 
3.16) package. 

The GEANT results were checked by comparison with the data in several 
ways. The total number of produced (5-electrons in "beam" (no interaction 
in the target) events and their pseudorapidity distribution were compared 
directly with the data for different targets and target thicknesses, and good 
agreement was found. For data with interactions in the target we compare the 
total number of 5-electrons seen by the Multiplicity Detector. We obtain this 
value from the data by extrapolation of the number of hits in the Multiplicity 
Detector Nl^^{Et) as a function of Et to the point -E* = (Fig. |) from the 
region Ef > 4 GeV, where N^j^^ grows linearly with Ef and the number of 
produced ^-electrons does not depend on Ef. The extrapolation gives for the 
number of 5-electrons values of about 13 and 24 for the 1% and 2% Au targets 
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respectively, which agree well with the GEANT simulation. The number 
of produced 5-electrons does not depend on Et in the region Et > 4 GeV 
because ^-electrons are mostly produced by the projectile nuclei before the 
collision due to the relatively large value of projectile charge Z, whereas the 
production of 5-electrons by final state particles is negligible. For very low 
Et events one expects an increase of (5-electron contribution to N^l^^ due to 
projectile fragmentation into high Z fragments, as can be seen in Fig. ^. The 
value of NJ^h^ at Et = corresponds to the number of 5-electrons produced 
by the projectile traversing the entire target. 

The GEANT (5-electron simulations have no free parameters which affect 
the results. The cutoff parameters in GEANT affect the low energy part of 
the (5-electron energy spectrum, but because of the absorption of low energy 
electrons in the aluminum absorber plates, we are insensitive to the choice 
of these parameters. The pseudorapidity distributions of (5-electrons used for 
the corrections for 1% Au target data are shown in Fig. ^ (curves 2 and 3). At 
the level of accuracy required for this analysis, he distribution of 5-electrons 
does not depend on the centrality of the collision. 

Gamma conversions. Another important correction is due to 7-conversions 
in the target and absorber. This correction was calculated using FRITIOF |Tl]| 
and RQMD |0 generated events combined with a GEANT simulation of the 
detector. In particular the ratio of pseudorapidity distributions of hits in the 
multiplicity detector due to 7-conversions and the pseudorapidity distribu- 
tion of charged particles was calculated for different centralities and target 
thicknesses. The results do not depend on the particular event generator 
used. They show that for a 1% Au target approximately 5% of the charged 
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multiplicity seen in the multiplicity detector is due to 7-conversions in the 
target and an additional 1% is due to conversions in the absorber. The effect 
depends slightly on the pseudorapidity: it is about 8% in the low part of the 
pseudorapidity region and about 4% for high pseudorapidities. About 30% 
of all e"'"e~-pairs occupy the same pad of the Multiplicity Detector. This 
effect results in the distortion of the Poisson statistics for the hit multiplic- 
ity distribution. The distortion caused by the two-particle correlations for 
produced secondaries is negligible, due to the very small magnitude of the 
correlations . The method used to correct for the 7-conversions is to apply 
the Poisson correction to the value calculated from the mean pad occupancy 
and then to subtract the distribution of electrons and positrons considering 
the pairs occupying the same pad as one charged particle. For the details of 
this procedure see Appendix A. The distribution used for the correction of 
the 1% Au target data is presented in Fig. ^ (curve 1) and has been used 
to obtain the charged particle multiplicity distribution for the highest cen- 
trality bin (see Fig. |^). The correction for other values of centrality scales 
approximately as dNch/dr]. 

4 Results 

4.1 Correlations among global variables 

The charged particle multiplicity measured in the Multiplicity Detector is 
strongly correlated with transverse energy deposited in TCal and/or PCal, 
and anti-correlated to the ("zero-degree") energy deposited in UCal. We 
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present these correlations in Fig. 0, where we plot the uncorrected (raw) 
number of hits in the Multiplicity Detector versus energy in each of the 
detectors. The relative abundance of the events with high multiplicity is an 
effect of trigger thresholds. The shapes of the the distributions in TCal/PCal 
Et for the events with fixed multiplicity, and the distributions in multiplicity 
for the events with fixed TCal/PCal Et, are close to Gaussian distributions. 
The length of the crosses in Fig. |^ indicate the widths of each distribution 
in different regions of the plot. 

The transverse energy deposited in the target (and/or central) region is 
strongly anticorrelated with the impact parameter of the collision. Therefore 
one can infer from Fig. that there are no drastic changes in the charged 
multiplicity fluctuations between central and non-central events. The width 
of the correlation between charged particle multiplicity and the PCal Et, 
which corresponds to the transverse energy deposited by all (charged and 
neutral) particles in the same pseudorapidity region, indicates that there 
are no large scale fluctuations in the ratio of charged and neutral particle 
multiplicities. 

For the current analysis the TCal Et data were used primarily as the mea- 
sure of centrality of the collisions. Note that these data are almost free from 
contamination by upstream interactions (except for low centrality events, 
where the correction for upstream interactions is small but not negligible). 
To compare our present multiplicity data with the published dEt/drj mea- 
surements in Au+Au collisions |Q, we use as the measure of centrality PCal 
Et data, which were effectively corrected for the leakage in the calorimeter 
through normalization to our previous measurements 0]. 
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The differential cross sections in TCal Et for different targets are shown 
in Fig. ^. The centrahty of the colhsions can be inferred using the ratio 
crtop{Et)/(Tgeom (shown on bottom plate of Fig. H), where atop{Et) is defined 
as 



and the geometrical cross section for the collision of A and B nuclei is (Jgeom = 
7r(i?A + RbY; Ra,b = 1.2A^/^ fm. The U, Cu, and Al targets data were 
obtained by (multilevel) triggering on PCal Ef. The width of the correlation 
between PCal Et and TCal Et is relatively large in comparison with the 
difference between different trigger thresholds. A special weighting technique 
was used to measure da/dEt{TCal) from data triggered with PCal Et at 
several thresholds (see Appendix B). 

The differential cross sections da/dEt (TCal, PCal) for the Au target 
agree well with our previous measurements [^]. Central collisions of Au 
nuclei with a U target produce about 20% more transverse energy in the 
TCal region than collisions with a Au target. However in the more central 
(PCal) region the produced transverse energy differs only by about 5%. We 
observe a similar difference between the two targets in the peak values of 
charged particle densities (see below). 

4.2 Charged particle pseudorapidity distributions 

Charged particle pseudorapidity distributions for different targets are pre- 
sented in Fig. ^ for different centralities. The centrality for the different 
TCal Et regions can be estimated from the ratio crtop{Et) / ageom shown in 




(8) 
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Fig. ^. For the discussion below it is important to note that for values of 
Et close to its maximum value the mean multiplicity depends weakly on the 
Et cut (for Au+Au collisions this is clearly evident in Figures ^ and |^ for 
TCal Et greater than 21 GeV). This suggests a simple criterion for defining 
the similar centralities for collisions between different projectile and target 
nuclei. The highest Et bins shown in Figures |^(a)-(c) are chosen to sat- 
isfy the requirements of a weak dependence of the multiplicity on the Et cut. 
These cuts correspond to values of <Jtop{Et) / (Tgeom from approximately 2% for 
Au+Al to 0.2% for Au+U, which correspond closely to cuts used to define 
central events in our studies of transverse energy production 

We have studied the effects of systematic errors in the charged particle 
pseudorapidity distributions as a function of pseudorapidity. From uncer- 
tainties in the corrections to the data (beam position, (5-rays, 7-conversions) 
we estimate the systematic error in the magnitude to be about 3% in the 
mid-pseudorapidity region for central Au+Au collisions. The uncertainty is 
slightly larger in the low pseudorapidity region (about 5%) and for lower 
centralities. We have also considered possible errors in the position of the 
peak found from Gaussian fits to the pseudorapidity distributions. These 
errors were evaluated by a variation of parameters in the distributions used 
for the corrections. It was found that the centroid of Gaussian fit is surpris- 
ingly stable for such variations. The systematic uncertainty was estimated 
to be not more than 0.03 units of pseudorapidity; it is mainly due to the un- 
certainty in the vertical position of the beam (which gives a contribution of 
about 0.02 units of pseudorapidity) and to the uncertainty in the correction 
for 7-conversions (also about 0.02 units). The uncertainties due to statistical 
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errors, are much smaller, than the systematic uncertainties. 

For central and mid-central collisions and heavy targets (Au and U) the 
pseudorapidity distributions are well fitted by Gaussians. We show the fits 
for the highest Ef regions in Fig. |^ (a) and Fig. ^ (b). If we compare the 
fits for both targets, we see that for the U target the height is about 6% 
larger (about 285 and 268, respectively). The position of the peak is shifted 
to a lower value of pseudorapidity (1.71 for U, and 1.76 for the Au target), 
consistent with naive expectations for the heavier target. The distribution 
has the same width (about 1.05) as for the Au target. Note that the widths 
of analogous distributions for very central collisions of a Si beam with Al, 
Cu, and Pb targets [@] have very nearly the same value. For gold collisions 
with Cu and Al targets (Fig. |^ (c)) the shape of the distributions becomes 
non-Gaussian. The peaks of the distributions are shifted to the higher values 
of pseudorapidity, as expected for collisions with light targets. 

For low centrality Et regions the pseudorapidity distributions for Au+Au 
and Au+U collisions exhibit an enhancement for low pseudorapidities. Monte 
Carlo studies imply that this asymmetry is caused by slow protons from 
target fragmentation. This hypothesis was checked by studying separately 
pseudorapidity distributions of hits with pulse heights around the minimum 
ionizing peak (where one does not expect contributions from slow protons). 
The distributions of low pulse height hits do not show such an asymmetry; 
thus the low pseudorapidity enhancement appears to be entirely due to hits 
with high pulse heights. 

We compare our results with the RQMD [g and FRITIOF ]ll|] event 



generator predictions in Fig. |T0| for the most central Au+Au collisions. The 
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centrality for Monte Carlo events was defined in accordance with the top 
cross section calculated using the transverse energy deposited in the TCal or 
PCal pseudorapidity region, as is done for the data. When using TCal Et 
as a measure of centrality both event generators underpredict the peak value 
of the pseudorapidity density (Fig. |10| (a)). FRITIOF also overpredicts the 
position of the peak 2.2 in comparison to the experimental value of 1.76). 
Similar trends were seen for FRITIOF results on pseudorapidity distributions 
of transverse energy [Q. In Fig. [lO| (b) we compare the event generator pre- 
dictions with the data using PCal Et as a measure of centrality. Whereas the 
difference between FRITIOF results and the data remains almost the same 
as in Fig. |l^ (a), the agreement between RQMD and data is significantly 
better. The origin of the different behavior is in the difference in the corre- 
lations between Nch with TCal Et and PCal Et for the data and the event 
generators. The FRITIOF Nch ^ TCal Et and Nch ^ PCal Et correlations are 
close to the observed ones, and the results of comparison of pseudorapidity 
distributions is insensitive to the choice of TCal or PCal Et as a measure of 
centrality. On the other hand RQMD exhibits a very tight Nch ^ PCal Et 
correlation, while the Nch ^ TCal Et correlation is rather loose. This results 
in lower mean Nch for the events with highest TCal Et in comparison with 
the mean value of Nch for events with highest PCal Et. Thus, for the same 
value of top cross section, the charged particle pseudorapidity densities for 
RQMD are quite different for the two cases. 
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4.3 Comparison with dEt/drj and evaluation of Et per 
charged particle 

Below we combine multiplicity data with our earlier measurements |Q of 
dEt/drj in Au+Au collisions at a similar, but not identical, energy. Due 
to the strong correlation between impact parameter and transverse energy 
produced in the central region (in our case PCal Et) it was found useful |^ to 
introduce the value the Et for an average collision with impact parameter 
h < 0.5 fm {o'top{E^)/(Tgeom = 0.22%), and use the ratio Et/ E^ as a measure 
of centrality. For our case this yields E^ 318 GeV. In Fig. 0(a) we present 
the charged particle pseudorapidity distributions for different PCal Et regions 
centered approximately at Et/E^ = 0.2, 0.4, 0.6, 0.8, and 1.0. The smooth 
curves represent Gaussian fits to the data. The peak position (r/peafc) and the 
width ((Teta) of Gaussiau fits are presented in Fig. |TT](b). The open circles 
in Fig. |ll|(b) show the dependence of parameters of Gaussian fit to dEt/drj 
distributions from Ref. 0]. The peak position of dNch/drj is very close to 
that for dEt/drj. It is smaller for lower centrality, presumably due to the 
relatively larger contribution of slow protons from the target. The width of 
the dNch/drj distribution decreases as centrality increases (similar to dN^h/dr] 
data for both the Au and Si beams), and is larger than the corresponding 
width of the dEt/drj distribution. 

We use the dEt/drj distributions for the same centrality (value of Et/E^) 
as shown in Fig. |Tl| to calculate Et per charged particle as a function of pseu- 
dorapidity and centrality, and to compare them with RQMD and FRITIOF 
predictions (see Fig. |1^). By Et per charged particle we mean the ratio of 
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average value of total transverse energy to the average charged particle multi- 
plicity Ef/Nch (which is not the mean transverse energy of charged particles). 
Note that the value of transverse energy per charged particle is less sensi- 
tive than the absolute pseudorapidity spectra to uncertainties in defining 
the centrality cuts for experimental and Monte Carlo generated events. The 
pseudorapidity dependence of Et per charged particle shown in Fig. |l2|(a) was 
calculated as a ratio of Gaussian fits to dEt/drj ^ and dNch/drj distributions. 
The observed large value of Et per charged particle about 0.75 GeV (for cen- 
tral collisions and in the central pseudorapidity region) is rather remarkable. 
If we take into account that the dEf/df] spectra were measured at slightly 
higher energy (11.4 GeV/nucleon compared to 10.8 GeV/nucleon for current 
data) and re-scale Et per charged particle with the available energy [1] , which 
differs for both cases by approximately 4%, we get the value of 0.72 GeV. 
This value is significantly higher than that in p+Pb (about 0.45 GeV |T3|), 
Si+Al and Si+Pb colhsions (0.55-0.59 GeV and 0.52-0.54 GeV, respectively 
[0,0) at an even higher beam energy of :^14.6 GeV/nucleon. Note that we 
compare the Et per charged particle for the central events triggered on Et. 



The pseudorapidity distributions in References and |T0[ were studied as a 
function of total charged multiplicity and cannot be used directly; but using 
the correlation between multiplicity and PCal Et it is possible to estimate 
the peak value of the distribution for collisions with the highest PCal Et. For 
example, for Si+Pb collisions this estimate gives a value of approximately 
115-120 particles per unit of pseudorapidity for central Si+Pb collisions. 
Combining this value with the peak value of dEx/dr] ^ 62 GeV observed in 
[^] one finds the value of Et per charged particle cited above. 
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The comparison of model calculations with the data shows that the 
RQMD event generator describes the pseudorapidity and nontrivial central- 
ity dependence of Ef per charged particle rather well, although the generator 
predicts the peak position to be at a slightly higher value of pseudorapid- 
ity. It cannot be excluded that this apparent disagreement is in part due 
to the representation of dNch/df] and dEt/drj as Gaussian functions and to 
the fact that the two data sets were taken at shghtly different beam ener- 
gies. FRITIOF does not reproduce the pseudorapidity dependence and, in 
disagreement with the data, shows no dependence on centrality. 

The possible origin of rather high value of Ef per charged particle was 
studied using the RQMD event generator, since it exhibits a dependence on 
centrality and pseudorapidity which is quite similar to the data . For this 
purpose Si-|-Pb at 14.6 GeV/c and Au-|-Au at 11.4 GeV/c collisions were 
studied. It was found that centrality dependence is almost totally due to 
changes in transverse energy deposited by nucleons. For less central events 
(and the collisions of nuclei of very different sizes, such as Si and Pb) the 
relative contributions of target spectators is rather large, which causes a 
decrease in Et per charged particle. The second reason for the difference in Et 
per charged particle between central collisions of light and heavy projectiles 
is a difference in the ratio of charged to all final state particles (caused by 
different relative numbers of protons and neutrons) . It was observed in this 
model, that the transverse energy of produced particles (pions, kaons) do not 
exhibit any strong dependence on incident energy (within the energy range 
considered), size of the target or projectile, and centrality of the collision. 
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5 Conclusion 



We have presented an analysis of charged particle multiplicity distributions 
in collisions of Au projectiles with Au, U, Cu, and Al targets with different 
centralities. The results are corrected for beam movement, upstream inter- 
actions, 5-ray production, and 7-conversions. In studying the correlation 
between charged multiplicity and energy deposited in the calorimeters, we 
observe little or no change in the fluctuations of charged multiplicity as a 
function of centrality and no large scale fluctuations between charged mul- 
tiplicity and transverse energy deposited in the same pseudorapidity region. 
The maximum value of the charged particle pseudorapidity density for very 
central Au-|-Au coUisions is close to 270 and about 5% larger for the U target. 
The transverse energy per charged particle grows with increasing centrality. 
For central Au-|-Au collisions it is close to 0.72 GeV, significantly higher than 
in p-|-Au or Si-|-Pb coUisions. 

The FRITIOF event generator underpredicts the number of produced 
particles in the central region. The peak of the FRITIOF pseudorapidity 
distribution occurs at a much larger value of pseudorapidity than in the data. 
RQMD also underpredicts the charged particle density if one selects central 
events using transverse energy in the target region. If transverse energy in 
the central pseudorapidity region is chosen as a measure of the centrality, 
the description is better. The reason for this hes in the looser correlation 
between multiplicity and transverse energy in the target fragmentation region 
in the events generated by RQMD, in comparison with the data. The RQMD 
description of the centrahty dependence of Et per charged particle (calculated 
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in the region 1.5 < f] < 2.0) is rather good; considering the pseudorapidity 
dependence of Ef per charged particle calculated for central collisions RQMD 
predicts the peak position of the distribution at slightly higher pseudorapidity 
value than data does. It is hkely that the degree of re-scattering is even more 
important in the multiparticlc production in heavy nucleus collision than it 
is implemented in the event generators; thus the production of thermalized 
hadron matter is more probable in such collisions. 
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Appendix A: Correction for 7-conversion 

Let us denote by n the mean charged particle multiphcity in a particular 
pad; n is the quantity of interest. The 7-conversions to e"'"e~-pairs result in 
additional hits in the pad characterized by the mean multiplicity of uncor- 
related (from different pairs) particles riji and mean number of pairs n^2, 
when both the e"*" and e~ from the pair occupy the same pad. In terms of 
these quantities the probability of a pad not being occupied (the input for 
the current analysis) is: 

po = 1 - e""-"^^-''^^ (9) 

The corresponding (to Poisson statistics) effective mean multiplicity equals: 

h ^ —\n{l — po) ^ n + n^i + ny2- (10) 

To obtain the value n one should subtract from n not the true mean number 
of e"*" and e~ from 7-conversions {ue = n^i + 2n-y2), but the number of 
"apparent" particles, considering the pairs occupying the same pad as one 
particle (n^i + 71^2) ■ 

The distortion to Poisson statistics for multiple hits in the same pad due 
to 7-conversions can be evaluated from the probabilities of single and double 
hits: 

pi = (n + n^i)e-"-'^^^(l-e-"-^), (11) 
P2 = {n + n^i)'e-"-"-^(l - e-"^2)/2 + n^2e-"-='(l - e""-"^!) (12) 
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Appendix B: Reconstruction of distribution from 
multilevel triggered data 



Below we discuss a technique which permits the construction of unbiased 
experimental distributions from multilevel triggered data for the case where 
the variable of interest is weakly correlated with the variable on which the 
data are triggered. We make use of the expression which gives the distribution 
dw /dx m. a, quantity integral over the distribution dw / dy in another 

quantity y, on which the trigger decision is based: 

dw f dw . . dP , . , 

where dP/dx{x;y) is the distribution in x for events at a fixed value of y. 
It is important to recognize that, whereas dw/dy is the distribution of the 
triggered quantity which is strongly influenced by downscaling (i.e. only 
a known fraction of the events above the trigger threshold are recorded), 
the function dP/dx{x]y) expresses the natural correlation between the two 
quantities and is independent of the trigger. If downscaling is introduced in 
the trigger, the distribution in y of the events which are written to tape is 
given by: 

dw dw , . 

where W{y) is the probability for the event with certain value of y to be 
written onto the tape and is constant within each trigger level region. We 
combine Eqs. (|lD-(0) to obtain 

S = /l4)^^'^^^"''^'^'' ^''^ 
25 



which can be used directly for the calculation of distributions of quantities 
from multilevel triggered data. Eq. ([T5|) has a very simple interpretation: 
in the calculation of distributions of some quantity x each event should be 
weighted with the inverse probability of the event being triggered. In other 
words, to calculate the natural distribution of events in x, we evaluate 

eventtTnAxW^d/)' ^ ^ 

where the sum is taken over all triggered events written on tape with value 
of X within the region Ax. 
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Figure Captions 



1. Experimental setup of E877 at BNL. For this analysis, we use data 
from the Multiplicity Detector, Target, and Participant calorimeters 
(see insert). 

2. E877 Multiplicity Detector, consisting of two identical 300 /xm thick 
silicon disks, segmented into 512 pads. The combined pseudorapidity 
coverage is 0.87 < rj < 2.46. 

3. Fits of the pulse height distributions by modified Landau distributions. 
Contribution from 1, 2, and 3 m.i.p. are shown. The dashed line 
indicates the Gaussian fit to the electronic noise. 

4. Distribution of the horizontal (solid histogram) and vertical (dashed 
histogram) coordinates of the hit centroid as a function of the beam 
position measured in the Beam Vertex Detector. Straight lines repre- 
sent linear fits to the distributions. 

5. The dependence of mean number of hits in the Multiplicity Detector 
as a function of TCal Et is shown for 1% (solid histogram) and 2% 
(dashed histogram) Au targets. The extrapolations to Et — gives an 
estimate of the numbers of 5-electrons. 

6. Charged particle pseudorapidity distributions used as corrections for 
Au+Au collisions (1% target). (1) Distribution due to 7-conversions 
for TCal Ef region 22-28 GeV, (2) distribution of (5-rays produced in 
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the target and in the absorbers, and (3)Distribution of 5-rays from 
upstream materiaL 

7. Transverse energy in TCal, PCal, and downstream energy deposited in 
UCal plotted versus the raw number of hits in the Multiphcity Detector. 
The crosses indicate the width (a) of the distributions in each variable. 

8. (Top panel): Differential cross sections in TCal Et for collisions of a Au 
projectile with U, Au, Cu, and Al targets. For clarity Au+Au, Au+Cu, 
and Au+Al cross sections are multiplied by factors of 0.1, 0.01, and 
0.001 respectively. (Bottom panel): The cross section integrated above 
Ef (see Eq. ^ and text thereafter), normalized to the geometrical cross 
section. 

9. Charged particle pseudorapidity distributions for different centralities 
(TCal Et regions). For the most central region we show also the Gaus- 
sian fit to the distribution, (a) Au+Au collisions; (b) Au+U collisions; 
(c) Au+Cu and Au+Al collisions. 

10. (1) Experimental charged particle pseudorapidity distribution for Au+Au 
collisions, compared with RQMD (2) and FRITIOF (3) predictions for 
the equivalent centrality; (a) TCal Et energy region 22-28 GeV ; (b) 
PCal Et energy region 310-330 GeV. 

11. (a): Charged particle pseudorapidity distributions in Au+Au collisions 
for the different centralities (PCal Et regions). The curves are Gaussian 
fits to the data, (b): Peak position and the width of the Gaussian fit 
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to dNch/drj as a function of centrality {Et/E^, for details, see text) in 
comparison with parameters of the analogous fits to dEt/drj 0], shown 
as open circles. 

12. Transverse energy per charged particle for Au+Au collisions, compared 
with predictions from RQMD and FRITIOF. The data are derived 
from Gaussian fits to dEt/drj taken at 11.4 GeV/nucleon and dNch/drj 
taken at 10.8 GeV/nucleon, as discussed in the text, (a) pseudorapidity 
dependence for central events (PCal Et energy region 310-330 GeV); 
(b) PCal Et dependence of Et per charged particle in the region 1.5 < 
rj < 2.0. 
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http://arXiv.org/ps/nucl-ex/9412003v2 



This figure "figl-2.png" is available in "png" format from: 



http://arXiv.org/ps/nucl-ex/9412003v2 



This figure "figl-3.png" is available in "png" format from: 



http://arXiv.org/ps/nucl-ex/9412003v2 



This figure "figl-4.png" is available in "png" format from: 



http://arXiv.org/ps/nucl-ex/9412003v2 



This figure "figl-5.png" is available in "png" format from: 



http://arXiv.org/ps/nucl-ex/9412003v2 



This figure "figl-6.png" is available in "png" format from: 



http://arXiv.org/ps/nucl-ex/9412003v2 



This figure "figl-7.png" is available in "png" format from: 



http://arXiv.org/ps/nucl-ex/9412003v2 



This figure "figl-8.png" is available in "png" format from: 



http://arXiv.org/ps/nucl-ex/9412003v2 



1 Si Multiplicity 
1 Detector \ 

! SI S2 S4 S3 \/ 1\ / 






1" ZJ-k 

BVer 1 BVer 2 TCal ^ 

' Si Detector / 
1 PCal 







Target 



TOFU 



Beam 




Target 

Scint TCal Magnet j^^^, 
units 



X 



Scint 



Forward 
Scint 



UCal 



34 
I 



meters 




Figure 2: 



ring 7, pod 33 




ADC channels 

Figure 3: 



10 

'o 




-7 -6 -5 -4 -3 -2 -1 



Xbver (mm) 

Figure 4: 




5 10 15 20 

Et (GeV) 




Figure 6: 



> 

20 

O 



UJ 



10 



4) 





400 



300 



200 



100 





2000 

CD 

LtJ 



1000 




50 100 150 200 250 300 350 

number of hits 



Figure 7: 



> 10^^ 



E 10 



X) 



1 ^ 
-1 = 

10 

-2 

10 k 

1 



e -1 
8 10 



10 



CL 

o 



b 



10 



X 1.0 



X 0.1 



xO.OI ^ 



X 0.001 



\ uH 

Au ^ 



Cu 

11 1 
1 1 til 








Al\ ICu Au^ \U 



10 15 20 25 30 

Et (GeV 



Figure 8: 



300 



ID 



250 



200 



150 



100 



50 



: Au+Au 






TCal E, (GeV) 






[22,28] 






[18,20] 






[13,15] 






[8,10] 






[4,6] 








300 



250 



200 



150 



100 



50 



1.5 



2.5 



7] 



Figure 9 (a) 




Au + U 



TCal E, (GeV) 

[25.5,28] 

[19,21] 

[14,16] 

[9,12] 
[4,7] 



1.5 



2.5 



V 



Figure 9 (b): 




Figure 9 (c): 




V 

(a) (b) 
Figure 10: 



^ 1 
> 

O 0. 9 

_q; 0.8 

o 

•-^0.7 
o 

CLO. 6 

"D 

^0.5 



CD 

Q- 



0.4 
0.3 ^ 
0.2 
0.1 






300 



250 - 



200 



150 - 



100 



(a) 



Figure 11: 









Au+Au 


L DATA 






--^^c;^. RQMD 




FRITIOF 







§0 

a; 
o 

= 

D 

Q-0 

|o 

o 



0) 

Q_ 



2 
(o) 





0.5 






(b) 






• -i>- 










'- -0- 




Au+Au 






• - DATA 






- RQMD 






□ - FRITIOF 









100 



V 



200 

(b) 



300 400 
Et (GeV) 



Figure 12: 



